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\ Abstract 

<D ■ 

Ph ■ We study a 1-D array of Josephson coupled superconducting grains with ki- 

netic inductance which dominates over the Josephson inductance. In this limit 
the dynamics of excess Cooper pairs in the array is described in terms of charge 
solitons, created by polarization of the grains. We analyze the dynamics of these 
topological excitations, which are dual to the fluxons in a long Josephson junction, 
\^ • using the continuum sine-Gordon model. We find that their classical relativistic 

motion leads to saturation branches in the I-V characteristic of the array. We then 
discuss the semi-classical quantization of the charge soliton, and show that it is 
consistent with the large kinetic inductance of the array. We study the dynam- 
ics of a quantum charge soliton in a ring-shaped array biased by an external flux 
through its center. If the dephasing length of the quantum charge soliton is larger 
than the circumference of the array, quantum phenomena like persistent current 
and coherent current oscillations are expected. As the characteristic width of the 
charge soliton is of the order of 100 /xm, it is a macroscopic quantum object. We 
discuss the dephasing mechanisms which can suppress the quantum behaviour of 
the charge soliton. 
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1 Introduction 

Arrays of Josephson junctions in 1-D, 2-D or 3-D have been studied exten- 
sively in recent years, both theoretically and experimentally [1J. When the 
capacitance of the junctions is small, the arrays are usually characterised by 
the Josephson energy, J2i Ej[l — cos(0j — 0j+i)], and by the charging energy, 
\ QiC^Qj. Here 0j and Qi denote the phase and the charge on the i'th 
grain of the array, respectively, is the inverse capacitance matrix, and 
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Figure 1: An equivalent electric circuit of a 1-D array of serially coupled 
Josephson junctions. 



Ej is the Josephson coupling energy. This description in terms of variables 
defined on the grains and not on the junctions is consistent with the fact that 
the kinetic and the geometric inductances of the grains are typically smaller 
than the Josephson inductance. As a result, the charge redistribution time 
in the grains is shorter than the tunneling time. In this paper we study the 
opposite limit, which is also experimentally accessible, namely a 1-D array 
where the kinetic inductance of the grains 
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dominates over the Josephson inductance 
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Here m* and e* are the Cooper pair mass and charge, n s the Cooper pairs 
density, l x the length of a grain and S the cross section of a grain. The large 
kinetic inductance means that in this case the charge redistribution time 
in the grains is longer than the tunneling time, thus the dynamic variables 
should be defined on the junctions of the array and not on the grains. This 
array can be represented by the electric circuit shown in Fig. [D. C denotes 
the self-capacitance of the superconducting grains, while the combined effect 
of the Josephson and charging energies of the junctions results in a non-linear 
capacitance, C, as we explain in the next section. We show that in this kind 
of array the concept of 'charge soliton' J2|-[]8|] arises, i.e. an excess Cooper 
pair in the array gives rise to a compact topological solitonic excitation. This 
appears to be in contrast to the usual model which does not incorporate the 
inductive effects. That model suggests that an excess Cooper pair delocalises 
as a consequence of the Josephson tunneling. We show, however, that a suf- 
ficiently large kinetic inductance decouples the individual junctions quantum 
mechanically. We study the dynamics of the charge soliton both classically 
and quantum mechanically. 

The paper is organised as follows: In section (2) we develop a continuum 
approximation of a serially coupled array of Josephson junctions with a dom- 
inant kinetic inductance. In section (3) we show that this array has compact 
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solitonic excitations ('charge solitons'), and discuss some of their classical 
properties and dynamics. In this section we discuss the small amplitude 
oscillations of the array ('plasmons') as well. In section (4) we study the 
classical dynamics of the charge soliton further, using collective coordinates. 
The quantization of the charge soliton is done in section (5). We discuss 
the meaning of the semi-classical quantization of the soliton, and study its 
quantum dynamics in a ring-shaped array. We demonstrate that quantum 
charge solitons can, in principle, exhibit quantum phenomena without clas- 
sical analogues, like persistent motion in response to an external flux and 
coherent current oscillations. We then discuss possible dephasing mecha- 
nisms of charge solitons, and address the effects caused by the discreteness 
of the array. We summarize our results in the concluding section (6). 

2 Kinetic Inductance Dominated 1-D Array 
of Serially Coupled Josephson Junctions 

2.1 The Lagrangian 

We consider a chain of iV identical superconducting grains (thus forming 
N—l Josephson junctions). The junctions are characterised by the Josephson 
coupling energy and by the charging energy scale 

(3) 

° 2C K J 

We assume that C ~ 10 -15 Farad, and that Ej is of the same order as 
Ec- The grains are capacitively coupled to a conducting substrate with a 
capacitance Co <C C, which we assume to be Cq ~ 10 -17 Farad. The energy 
scale of this coupling energy, 

Ec = fi£ (4) 
is thus much larger than the junction charging energy 

E Co » E c . (5) 
The grains are characterised by the inductive energy scale associated with 

Lkin 

El = 7yT^ > (6) 

where $o = h/2e. As we have said in the introduction, we assume that the 
kinetic inductance of the grain dominates over the Josephson inductance. 
In fact, due to the numerical coefficient (2ir) 2 /2 difference in the relations 
(@) and (|), L kin should be larger than 2n 2 Lj for the inductive effects to 
be important. For a typical Ej of the order of 100 \ieV it means that L kin 
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dominates if it is 10~ 7 Henry or larger. This situation can be achieved, for 
instance, when l x w 10 fim and S ~ 10 3 nm 2 . Nevertheless we assume that 
the width of the grains is of the order of the London penetration depth to 
avoid tunneling of flux quanta through the grains. The width of the junctions, 
d, is much smaller than l x (typically d ~ 2nm), and the distance between 
adjacent grains (the unit cell) is denoted by a (a = l x + d). L = Na is the 
total length of the chain. We assume that the chain is very long (N 3> 1). 

Using the values given above, we find that the zero frequency impedance 
of a unit cell, 

ZlC = \J Lkin/ Co , (7) 

is of the order of 100 KQ, i.e. it is much larger than the quantum resistance, 
Rq = V(2e) 2 : 

Zlc > Rq ■ (8) 

Note that this impedance inequality can be expressed alternatively as an 
inequality of the coupling energy and the inductive energy scales 

E Co > E L . (9) 

A similar condition to (|]) has been studied before in the context of single 
electron tunneling in a normal junction |§, and it has been shown that it leads 
to a quantum mechanical decoupling of the junction from its environment. 
Using the same reasoning here, we are led to the conclusion that the condition 
(§) means that each junction is quantum mechanically decoupled from its 
environment, i.e. from the other junctions of the array. We can thus solve 
the Schrodinger equation for each junction separately, and obtain a local 
potential energy of the array. This situation has been named the 'local rule' 
in the context of single electron tunneling [|TUJ. 

The eigenstates of the junction i depend on q iy the dimensionless charge 
(in units of 2e) induced on this junction. As a function of q~i, the energy levels 
are made of a set of charging energy parabolas, with gaps at the intersection 



regions due to the Josephson energy [|TTJ-|jT4| (see Fig. |2]). The energy levels 



are, thus, periodic functions of with a period 1. Under appropriate condi- 
tions (not too small gaps, adiabatic changes) Zener transitions between the 
levels can be avoided fll5] , flEfl . We also ignore, for the time being, quasi- 



particle tunneling, which is a dissipative process. We discuss this issue in 
section (5). We thus may consider only the first level, which we denote by 
Eg.. This level represents coherent superposition of charge states in the bulk 
superconductors, which differ by one Cooper pair. Eg. is formally given as 
an eigenvalue of Mathieu's equation. As it does not have a simple analytical 
form when Ec is of the same order of Ej, and our results do not depend 
qualitatively on the exact form of E^, we adopt the following form 

Eg^-^Ecll-cos^Trq,)] . (10) 
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Figure 2: Energy levels of a Josephson junction as a function of $ for the 
case £7 j ~ i^c- 

The form fllOD preserves the correct parabolic dependence for small qi, and 
reduces the amplitude of the energy level from its maximal height (in the 
limiting case where Ej = 0) by a factor of 7r 2 /4. We emphasize that the 
important feature of Eg. is its periodicity, which allows us to represent the 
Josephson junction as a non-linear capacitor (see Fig. p]). In the next section 
we show that the periodicity gives rise to the soliton description. 

Due to the tunneling of Cooper pairs the variable qi is compact, i.e. 
qi + 1 = q~i. It is convenient to introduce an extended variable qi, which is 
the dimensionless charge (in units of 2e) brought to the i'th junction, is 
related to q~i through 

N 

qi = qi+ Qv ' ( n ) 

i'=i+l 

where Qi is the net charge on the i'th grain. Qi has, of course, only discrete 
values, while qi and q\ are continuous. This change of variables corresponds 
to changing from a 'reduced zone' scheme to an 'extended zone' scheme in 
the junction's energy bands (Fig. 0). This variable was used in the study 
of 1-D array of serially coupled normal junction as well @], ||]. In the next 
section we show the importance of qi for the solitonic description. The form 
of the energy of the junction (the potential energy) does not change when 
expressed as a function of qi 

E pot = E qi = j^E c [l - cos(2vrg 4 )] . (12) 
The voltage across the junction, V gi , is given by the derivative of the energy 
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levels with respect to the charge 



1 dE, 



V* = • (13) 



91 2e % 
Using ([12]) we express the voltage as 

V <H = ±V c sin2nq i , (14) 

where Vc = ff- 

Since is defined on the junction it is already contains an averaging over 
the fast tunneling process. A time dependent is therefore related to the 
slow process of charge redistribution in the grains by means of a supercurrent. 
This gives rise to an inductive energy in the grains, which serves as the kinetic 
energy of the array: 

E k in = -{2e) 2 L kin q} . (15) 

In the parameters range we consider, the kinetic energy scale is smaller than 
the potential energy one 

E C >E L . (16) 

The three inequalities: (|5|), (0) and ([H]) can be combined into a single con- 
dition for the energy scales of the system 

E Co >E C >E L . (17) 

The relation between the dynamic variable qi and the voltage Vi between 
the i'th grain and the substrate can be found by consecutive applications of 
Gauss' law 

ft = ?i - y c ° E Vi ' ' ( 18 ) 

i'=l 

where q± is the charge that was brought to the first junction of the array. 
From now on we assume that the continuum limit can be taken. (We will 
show the necessary condition for this soon.) Discreteness effects are discussed 
in section (5). In the continuum limit Eqs. ( p]) and (|18D have the form 

q ( x ) = q( x ) + £Q^)d^/a, (19) 

q(x) = q(0) ~ ^C f X V{Odi/a. (20) 
le Jo 

The array is thus described by the charge field q(x). The relation ( pUP be- 
tween q(x) and V(x) can be expressed in a local form 

V(x) = -a^q x (x) . (21) 
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We see that the q x (x) is the dimensionless charge between the grains and the 
substrate. The charging energy which couples the unit cells of the array can 
be expressed, therefore, as 

TP — n 2 ^jiLn 2 (99> 

-^coupling ~ « ^(J ^ ' — ' 



As we have mentioned above, its energy scale is E Co (|J). When Co < C 
we have Ec 3> Ec- In this case even small amounts of charge induce high 
voltages on the capacitors between the grains and the substrate, and these 
voltages strongly couple the Josephson junctions. In the opposite case when 
C is large, there is almost no voltage on the capacitors and the junctions are 
practically decoupled. A small Co is thus needed for the picture of serially 
coupled Josephson junctions. 

From the above discussion we conclude that the array we consider is char- 
acterised by the three energies: the potential energy (|i~2"|), the kinetic (or in- 
ductive) energy fll5|) and the coupling (or charging) energy (|22|) . When these 
three energies are combined, we get the following sine-Gordon Lagrangian 

C = i-(2e) W - ogrf - i^Ml,! - co S (2.,)] . (23) 

This is a novel description of a 1-D Josephson junctions array, which is valid 
when the condition ( |17"D holds. The three effects of the large kinetic induc- 
tance are reflected in the Lagrangian fl2"3"|): 1. An additional inductive energy, 
which is an inertial term. 2. A representation of each junction by a periodic 
charging energy, as a result of the quantum mechanical decoupling of the 
junctions. 3. A description of the array by degrees of freedom which are de- 
fined on the junctions and not on the grains. The Lagrangian ( P3"D is electro- 
magnetically dual to the Lagrangian representing a long Josephson junction. 
The latter case can be understood as the continuum version of an array of 
parallely coupled Josephson junctions. Interchanging parallel coupling with 
series coupling and inductors with capacitors one gets the Lagrangian of the 
serially coupled Josephson junctions. Note, especially that the periodic in- 
ductive energy in the long Josephson junction (i.e. the Josephson energy) is 
replaced here by the periodic charging energy. 



2.2 The Equation of Motion and the Hamiltonian 



Following the standard sine-Gordon treatment ]T7|], we redefine the 

charge field: q(x) — ■> q'(x) = q(x)/27i, and express the Lagrangian ([2^) as 

hv c 



C 



27T/3 2 



2vl 



Q 



— 2"(1 — cosg) 



c 



(24) 



The three bulk parameters: C, Lk in and Co are replaced in (^|) by Ac, vc 
and (3 2 . Here 

A c 



(25) 
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is the characteristic length of the system. The condition needed for the 
validity of the continuum limit is therefore 



Ac > a , (26) 

or 

C » Co , (27) 

which is consistent with the limit (|). This is another manifestation of what 
we have discussed above: a small Cq implies a large coupling, hence a large 
A c . Using the values given above we get Ac « 100 jum. The second param- 
eter in the Lagrangian (|24|), 

vc = , (28) 

is the wave velocity of the system. It is of the order of 10 _1 — 10~ 2 c, where c is 
the vacuum light velocity. It is related to Ac via the characteristic frequency 



° Ac \j L kin C ' ^ ^ 

which is of the order of 10 11 sec _1 . The third parameter in the Lagrangian 

(El), 

sets the energy scale of the system. It does not affect the classical equation of 
motion, but its value is important in determining whether the system behaves 
classical or quantum mechanically. We return to this point in section (5), 
where we discuss the quantum dynamics of the system. 

The equation of motion derived from the Lagrangian (p4|) is 
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— q - q xx + -rYsmq = . (31) 

v c A c 

It is a voltage equation for the junction, as can be shown more clearly by 
multiplying it by 2et> c Z/fci n /27r and using Eq. fl2"T|) to obtain 

1 1 2 2e 1 

—2eL kin q - —a —q xx + —V c smq = . (32) 

Z7T Z7T Go 27T 

The first term is an inductive voltage induced along the grains when the 
current is time-dependent. From Eq. (|2lD we see that the second term is the 
continuum form of Vi+i — Vi, i.e. it is the difference of the voltages between 
two adjacent cells and the substrate. The third term is the voltage across 
the junctions, resulting from the superposition of charge states (Eq. (fH|)). 
The voltage equation (|32|) is thus a Kirchoff's law for a closed loop of the 
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equivalent electrical circuit of the array (Fig. [I]). The conjugate momentum 
of the field q 

fir 1 (2e\ 2 
L ' ' ~ (33) 



J kinQ 



-= hn. 



dq a \2ix , 

is the number of flux quanta per unit length that has tunneled through the 
junctions of the array. Using n$ we get the Hamiltonian of of the system 



H 




2 qx+ A 2 C 



[1 — cos q) 



dx 



(34) 



When the array is coupled to an external voltage, V ext , the equation of 
motion (|3T| ) changes to 

1 1 



— q - q xx + jj sin q 



J c 



2n — — V ce u 
a z 2e 



where 



V C ell = — Vext 



(35) 



(36) 



is the part of the external voltage that is distributed on one unit cell. Equa- 
tion ([35|) represents, alternatively, the case where the array has a shape 
of a ring and an external flux is applied through its center. In this case 
V ex t = —&ext is the electromotiv force acting on the array. The flux source 
has, of course, the advantage that the effects of the leads are eliminated. In 
any case, equation (|35| ) can be derived from the following Hamiltonian 



H 




2n(3 -(n* -n 4e J + ^-^ 



:C + tt( 1 - cos ^) 



A 



c 



dx 



(37) 

In the case of a voltage source n<s> ext is defined as the integral of the external 
voltage per unit length and unit flux 

1 



Vext dt 



(38) 



while in the case of a flux source it is simply the dimensionless flux density. 
The external source thus appears in the Hamiltonian as a time-dependent 
gauge potential, in analogy to the external current in the long Josephson 
junction Hamiltonian [19]. The gauge nature of the external voltage gives 



rise to the following shift of the conjugate momentum 

hn<s> = ~ (tt) L kin q + Tin$ ext . 
a \2ty/ 



(39) 



Dissipation processes in the system produce additional q - dependent 
voltage drops. Ohmic dissipation can be represented phenomenologically by 
adding to each unit cell a resistor connected to the other elements in this cell 
in series. In this case the voltage equation (R5I) becomes 



— q + -^C Rq 



q x 



'c 



+ - T smq 



A 



c 



2>K-z — V ce ii 
a 2 2e 



(40) 
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This representation, which was named the 'serially resistive junction' (SRJ) 
in ||, is the analogue of the RSJ model |20fl , 

3 Charge Solitons and Plasmons 
3.1 A Static Charge Soliton 

Since the 1-D array of serially coupled Josephson junctions can be described 
by a sine-Gordon Lagrangian fl2"4]), we expect that it has solitonic excitations, 
i.e. compact, stable topological configurations. Using the definition of q as 
an extended variable (fL9f) , we observe that q(x) and q(x + 2ti) can be distin- 
guished if there is an excess or a deficiency of Cooper pairs in intermediate 
grains. The one soliton excitation represents the charging of the junctions 
(or the polarization of the grains) due to an excess Cooper pair in the array, 
and is called a 'charge soliton'. This term was coined in in the context of 
a 1-D array of normal tunnel junctions. Recently, charge solitons in a 1-D 
array of SQUID's were studied experimentally ||, and a zero current state 
below a threshold voltage was found. This voltage was interpreted as an 
injection voltage for a charge soliton. 

The charge soliton solution of Eq. (j31|) with the appropriate boundary 
conditions is (see Fig. ||) 

q so i{x) = 4 tan 

Its center is at X , which we take in this section to be zero. The excess 
charge of the Cooper pair is the topological charge of this soliton 

Q = 2eJ d x q sol dx = -2e . (42) 

We would like to emphasize once more that under the condition ([lj]) we con- 
sider here, the existence of a topological solitonic excitation and its stability 
do not depend on the exact form of the potential energy of the junctions, 
but only on its having degenerate minima. Thus our qualitative results are 
valid for other forms of the potential as well. 

As was mentioned above, charge solitons in 1-D arrays of normal tunnel 
junctions have been studied previously 0-(7||. In this context a question was 
raised whether a charge soliton can be regarded as a coherent dynamic object 
whose equation of motion contains an inertial term, as was proposed in Refs. 
@, || and HJ], or that it merely represents a static charge distribution profile, 
as was argued in Refs. []5| and ||. Here we have shown that this question 
should not rise in the Josephson junction array context. The coherence of the 
charge soliton ensues from the coherent superconducting ground state, and 
the inertia term comes from the kinetic inductance of the grains. Moreover, 
we have shown that the impedance condition @ should be met in order that 



exp 



X -A C 
Ar 



2tt 



(41) 
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Figure 3: The charge soliton configuration representing an excess Cooper 
pair in the array. The center of the soliton is taken to be X$ = 0. 



the concept of a charge soliton will be different from that of a point charge 
(be it a Cooper pair or an electron). 

From Eq. (f4"T|) we see that the characteristic length scale of the array, Ac, 
is the characteristic width of the soliton as well. In order to interpret the 
charge soliton as a particle its width should be much smaller than the total 
length of the array, i.e. 

L > A c • (43) 

This assumption is met when L > 10 3 fim. Here we assume that L « 10 3 \xm. 
The number of grains the soliton is spread over is 



N c = Ac/a = JC/C . (44) 



Nc is larger than one due to the continuum limit condition (|26|). For the 
parameters given above Nc = 10. when the condition (|26|) fails, one should 
take into account corrections to continuum sine-Gordon model. We address 
this point in section (5). The finite width of the charge soliton is clearly seen 
from its density, which according to Eq. (pll), is proportional to the profile 
of the voltage between the array and the substrate (see Fig. |j) 

v{x) = = ^iv^ sech {-B • (45) 

Ac sets the scale for the static distribution of voltages on the junctions of 
the array as well. Using Eqs. flUD and (31) we find that this distribution is 



proportional to the second derivative of the soliton configuration (see Fig. |5|) 
V q (x) = ^-V C A 2 C d xx q sol (x) = -—V c sech (-?-) tanh . (46) 
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Figure 4: The profile of the voltage between the array and the substrate 
induced by the charge soliton. V is measured in mV. 



The energy needed to create a charge soliton is the value of the Hamilto- 
nian calculated for a static solution (Eq. (f41"D) 



8 hv c 8 (2e) 2 16 

Eo - a~ C 2^ - wyvm ~ w? EcNc ■ (47) 

This rest energy depends on C and C but not on L kin , since it is determined 
by the potential and coupling energies. It can be written as the potential 
energy density (ec = Ec/S), times the effective area of the soliton (S e ff = 
SN C ) 

16 

Dividing Eq. ( (47| ) by v c we get the soliton rest mass 



E ~ 7nZ\2 e cSeff ■ 



Mo = E /v c = 7^(2e) 2 ^-^ . (49) 
(2ti) z a A c 

In analogy to the rest mass of a fluxon in a long Josephson junction [19], the 



charge soliton's rest mass is proportional to the inductance per unit length 
and inversely proportional to the characteristic length, Ac- Using the typical 
parameters we find that the charge soliton mass is of the order of 10 -36 Kg, 
i.e. six orders of magnitude less than the electron rest mass. This result 
indicates that the charge soliton should not be understood as a Cooper pair 
dressed with a polarization cloud, but as the polarization cloud itself. We 
return to this point when we discuss the dynamics of the charge soliton in 
the next section. 
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Figure 5: The distribution of voltages on the junctions of the array corre- 
sponding to a charge soliton configuration. V q is measured in mV. 



3.2 A Dynamic Soliton 

In order to describe a charge soliton moving with a velocity v, we make 
use of the Lorentz invariance of the Lagrangian to perform a Lorentz 



transformation of the static configuration (f|l]) and obtain 

'x — Xq — vt 



q so i{x,t) = q so i[l(x -vt)] = 4 tan <^ exp 



7 



A 



c 



(50) 



where 7 = 1/ ^Jl — (v 2 /vq). We thus expect that a relativistic charge soliton 
suffers a Lorentz contraction. Since the light velocity in the array, vc, is 
smaller than the vacuum light velocity, relativistic effects of the charge soliton 
can be observed more easily than relativistic effects of electrons or Cooper 
pairs. 

A moving charge soliton induces, of course, a current along the array. 
The spatial distribution of the current is given by 

2e 

H x ) = • ( 51 ) 

This is a current pulse with a width Ac, concentrated around the moving 
center of the charge soliton. It has the same profile as the voltage between 
the array and the substrate (see Fig. [|). The average current produced by 
the moving soliton is 

If 1 
I = - J I( x )dx = --2ev . (52) 

For a soliton moving with a velocity 10 6 m/sec, it is of the order of 0.1 nA. 
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3.3 Plasmons 



Besides topological solitons, the sine-Gordon Lagrangian (|24[) admits small 
amplitude excitations. Their dynamics is governed by the linearized equation 
of motion 

1 1 

—q ~ Qxx + -ryg = . (53) 
v c A c 

As this equation describes electromagnetic field oscillations with a confining 
potential, its solutions are longitudinal plasma oscillations ('plasmons') prop- 
agating along the array. The propagation of the plasmons does not involve 
any tunneling process. The plasmons have the dispersion relation 

to 2 = uj 2 c + v 2 c k 2 , (54) 

i.e. there is an energy gap hu>c in their spectrum with the corresponding 
temperature T g « 1 K. The plasmons have, therefore, a mass 

M P =*£= * (55) 
v 2 c A c v c 

which is of the order of 10 -37 Kg. The ratio between the mass of a plasmon 
to the mass of the soliton ( ^9|) is 2nj3 2 /8, i.e. it is of the order of (3 2 . 

Plasmons can also be excited when there is a soliton in the array. In that 
case they can be considered as vibrations of the soliton. Their analytical 
form can be found by expanding q around the soliton solution (|4T|) 



q(x, t) = q so i{x) + ipk{x) exp(iu k t) . (56) 

Substituting ( [56"D in the equation of motion ([H]), and linearizing with respect 
to if) one gets 

^k( x ) ~ (tanh- ikAc J exp(ikx) . (57) 

The dispersion relation is the same as above (Eq. (0)), but there exists 
now an additional zero mode (whose u = 0). It reflects the translational 
invariance of the system, i.e. the homogeneity of the array (at distances 
larger than a). As a consequence, the zero mode is proportional to the 
spatial derivative of the soliton configuration. 



4 Collective Coordinates for the Charge Soli- 
ton 

4.1 Definition and Equations of motion 

The topological stability of the charge soliton and its finite width allow for 
its interpretation as a particle. Thus we would like to describe the charge 
soliton by a pair of conjugate coordinates which correspond to its center of 
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mass, X, and to its momentum, P. This can be done by using the 'collective 
coordinates' method. This method was studied extensively in the context of 



general soliton theory [22]-[25], as well as for Josephson junctions in particular 
0i Hl> EH; H9' [p7|1 . The basic idea behind it is that the field q can be 
expanded around the soliton solution in a form similar to (Eq. ([56])) 

q(x, t) = q so i(x - X(t)) + c k i) k (x - X(t)) exp(iw k t) , (58) 

k 

where the c k are the normal modes around the soliton configuration, i.e. the 
plasmons, and the sum now does not include the zero mode. The soliton's 
center of mass is thus elevated into the role of a dynamic variable at the 
expense of the non-physical zero mode. The Hamiltonian ( |34| ) can be now 
canonically transformed into a form involving the new variables. Alterna- 
tively, in order to study the dynamics of the soliton itself, one can retain 
only the collective coordinate by assuming the form 

q(x,t) =q sol (x-X(t)) , (59) 

which means that the soliton is considered to be a rigid object moving with 
velocity X. This assumption is justified when the temperature is much lower 
than the plasmons' energy gap. The plasmons are then treated as a pertur- 
bation. In this method the particle interpretation of the soliton is clearly 
seen. 

The collective coordinates can be expressed in an explicit form by using 
the soliton density, d x q so i, as a weight function [p7f| . The center of mass of 
the charge soliton (the position of the excess Cooper pair) is given by 

X ^ ~~ 2tt / xdxqso1 dx ' ( 60 ) 
and the conjugate collective momentum by 

P^Jn q d x q sol dx. (61) 

One can check that X and P are indeed canonical variables by calculating 
their Poisson brackets: [X, P]pb = 1- Inserting the soliton configuration (|50| ) 
into the definitions ( |60|) and (|6lD we get the equations of motion of a free 
relativistic particle 

X = X + vt , (62) 
X = v , (63) 
P = lM X , (64) 
P = . (65) 
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4.2 The Dynamic Mass and the Hamiltonian 

The mass that appears in ( |64]) is actually the dynamic mass of the charge 
soliton 

'2e\ 2 f 

-) L kin J d x q sol dx . (66) 

Its value is identical to the rest mass (|49|) in the limit L 3> Ac, and differs 
from it by a factor of 2 in the opposite limit L <C Ac- As we consider here 
the first limit, we denote it by M as well. We can understand the origin 
of the dynamic mass by observing the way the charge soliton propagates. 
Starting from the static distribution of voltages on the junctions (Fig. |5|), 
the center of the charge soliton moves from its position in the middle of 
a grain towards one of the neighbouring junctions, say the right one, by a 
charge redistribution in the grains. A superposition of charge states in the 
two adjacent grains is built, and the (negative) voltage on this junction is 
reduced. When the superposition is of states of equal weight, the voltage is 
zero. As the motion continues, the charge redistribution increases the weight 
of the charge state on the right grain and the voltage on the junction is 
increased. When the absolute value of this voltage reaches the initial one, 
the center of the charge soliton has been shifted by one unit cell, i.e. it is in 
the middle of the right grain. One sees that the propagation of the charge 
soliton is determined by the kinetic inductance and not by the Josephson 
one. The dynamic mass leads us, therefore, to the same conclusion that 
we got from the rest mass: the charge soliton is the polarization cloud that 
accompanies the excess Cooper pair that exists in the array. 

Transforming now the Hamiltonian (34) into collective coordinates form, 
we get 



H = ^Mlv% + P*v c , (67) 
so the energy of the moving soliton is 

E = jM v c = 7^0 • (68) 

If we assume the non-relativistic limit, i.e. v <C v c , the Hamiltonian describ- 
ing the soliton particle reads 

H = M ° v ° + & • m 

where now 

P = M X . (70) 

The rest energy term in the Hamiltonian (^9|) is made out of the two charging 
energies (the last two terms in fl3"4])), while the contribution to the kinetic 



term in (j69|) , comes only from the inductive energy (the first term in (|34])). 
We thus see that the inductive energy, although being the smallest energy in 
the system, is the one that governs the dynamics of the charge soliton. The 
independence of the Hamiltonian (^) on X is another manifestation of the 
translation invariance of the system. 
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4.3 A Voltage Biased Array 

The collective coordinates can be used to describe a voltage (or a time- varying 
flux) biased array as well. Introducing the external voltage in the form 



&ext = —Vext , 

we find that the collective momentum is shifted to 



P = M X + 



ext 



L $ 

and the non-relativistic particle Hamiltonian is 

2irK& ext * 



(71) 



(72) 



H = M Vc + 



1 



P 



2M V L $ 
The equations of motion derived from (^) are 

2nh $ ext \ 



(73) 




L $ 



o 



and 



P = . 

Combining the two equations we get 

2ith $ ext 



M X 



L $ 



o 



(74) 



(75) 



(76) 



i.e. the external voltage accelerates the charge soliton. The origin of this 
acceleration is simply the electrostatic force exerted on the excess Cooper 
pair by the external voltage. In order that the rigid soliton assumption will 
be valid in this well, the external flux must be changed adiabatically, 

or the external voltage should be small enough 



$ext 




Vext 


$ 




$ 



(77) 



which means that V ext should be of the order of 10 fiV or less. 

When there are Ohmic dissipation processes in the array an application 
of an external voltage results in a steady state velocity (or current) of the 
soliton. The steady state is reached when the power gained by the voltage 
is equal to the power lost via the dissipation. Using the Hamiltonian (p4|). 
the equation of motion (0), and the average current ([5^), we find that the 
steady state condition is 



V P ~ 



ext 



Reffh 



(75 
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steady 



Figure 6: I-V characteristic of a voltage-biased dissipative array. Each branch 
corresponds to a certain number of charge solitons in the system. V ex t is 
measured in [iV and I s teady is measured in nA. The parameters are R = 10Q 
and vc = 10 _2 c. 

where the effective resistance of the array is constant in the non-relativistic 

case 

Reff = Wy^ R (79) 

and is / siea( iy-dependent in the relativistic case 



R 



■eff, rel 



\1 steady , 



L 1 



(27r) 2 aA c 



1 {I steady) R 



(80) 



The effective non-relativistic resistance of the array is thus increased by about 
two orders of magnitude, while relativity increases it further by the 7 factor. 
The I-V characteristic of the array is expected to show saturation branches, 
where each branch corresponds to a certain number of solitons reaching the 
limit velocity, vq (Fig. ^). 



5 Quantum Dynamics of the Charge Soliton 
5.1 The Semi- Classical Expansion 

In this section we study the quantum dynamics of the charge soliton as a par- 
ticle. For this we utilise the semi-classical quantization of the sine-Gordon 
theory |l7j , |L8| , [2^] , [^8| . The expansion parameter is the coupling constant 
/3 2 , which was defined in fl30|). In this method the total Fock space is taken 
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to be composed of disconnected sectors, each one corresponds to different 
topological boundary conditions, i.e. to a different number of solitons in 
the system. The ground state of each sector is the corresponding solitonic 
configuration. Here we concentrate on the one-soliton sector. Due to the 
translational invariance of the system there is, in fact, a degenerate family of 
eigenstates of the position operator, connected by space translations. Higher 
states are found by a semi-classical expansion around the ground state. The 
excitations of the first order correspond to the plasmons, and their quantum 
interpretation is as light particles scattering from the static massive soliton. 
These plasmons are, thus, the fundamental quanta of the theory. The degen- 
eracy of the states is completely removed in the second order, as the position 
eigenstates are replaced by momentum eigenstates, and the translation in- 
variance of the theory is recovered on the quantum level. The semi-classical 
expansion breaks down when (3 2 > 2, where the soliton becomes lighter than 
the plasmons. The soliton then takes the role of the fundamental quantum, 
and loses its correspondence to the classical particle configuration (in a sense 
it becomes 'too' quantum). Since the typical value of j3 2 is 10 _1 , we can use 
the expansion for the array. The parameter (5 2 can be expressed in the form 

= J* (81) 
Comparing Eq. flgTp with (H) and (0), we see that the condition for using 




the semi-classical expansion, (3 2 <C 1 is identical to the impedance condition. 
This is not a surprise, as the impedance condition is the one that enables us 
to decouple the junctions quantum mechanically. Our model of the charge 
soliton as a classical configuration is thus self-consistent. 

However there are several differences between the system we study and 
the field theoretical model. First of all, the array is very long (compare to 
Ac), but finite. Apart from a slight distortion to the soliton's shape that we 
neglect, the finiteness means that solitons can enter and leave the array, and 
also get reflected from the edges. To avoid this situation, we consider a ring- 
shaped array. Second, since the gap in the plasmons' spectrum is of the order 
of one Kelvin, their population can be made negligible if the temperature is 
kept below the gap. Thus we can discard all the plasmons' contribution to the 
dynamics. This assumption is equivalent to the rigid soliton assumption (|59D . 
A finite population of plasmons can be considered as an internal environment 



which produces a phase breaking mechanism [30|. We comment on this 
dephasing process at the end of this section. Another different feature, is 
that we couple the array to an external (classical) flux source as a gauge 
coupling, and study the quantum dynamics of the soliton in response to this 
source. Finally, the array we study deviates from the ideal sine-Gordon model 
by its discreteness, by the exact form of the potential energy, by structural 
inhomogeneities and disorder and by quasi-particle tunneling. The effects of 
these deviations from the ideal model are discussed below. 
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5.2 Persistent Motion of the Charge Soliton 

In the presence of an external flux, $> ex t, the assumption of rigidity leads to 
the following non-relativistic quantum Hamiltonian for a ring-shaped array 
of serially coupled Josephson junctions 

*=**+M p -™%)' ■ (82) 

Higher orders contributions to the energy give rise to quantum corrections 



to the soliton's rest mass |JjJ . The renormalized mass in the array language 
(up to the order of (3°) is 

M 0ren = ^-4r~ f 1 - ?1 = M o (l - ^) ■ (83) 



Ac 2n/3 2 v c 

However since ft 2 is small we can use M instead of M 0ren . As we have dis- 
cussed in the previous section, the Hamiltonian is X independent due to the 
homogeneity of the array. Thus it commutes with the collective momentum 
operator, P, and the eigenstates are collective momentum eigenstates with a 
discrete set of eigenvalues, pjq = hk N determined by the periodic boundary 
conditions 

k N = ^—N- N = 0, ±1, ±2, ... . (84) 

Lj 

The energy spectrum is discrete, too, and is given by (neglecting the constant 
term MqVq) 

e ^M&^ n -^ ? - (85) 

Defining an effective inductance by 

L ^ M °(S§) 2 (86) 

{L e ff « 10~ 5 Henry), the energy levels can be expressed in the form of 
inductive levels 

E N =-—{$ Q N-<S> ext ) 2 . (87) 

ZL eff 

The inductive form of the energy levels suggests the interpretation of as the 
number of flux quanta that has tunneled outside or inside the ring through 
one of the junctions. The quantization of P is, therefore, the statement that 
only an integral number of flux quanta can tunnel in or out of the ring. 
However, the conservation of the momentum means that there can be no 
flux tunneling in an homogeneous array, i.e. the external flux is completely 
screened. 

The spectrum of the charge soliton's Hamiltonian (]85|) is periodic with 
respect to the external flux with a period $o- It is composed of a set of 
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parabolas centered at $ eX f = iV$o- Each parabola intersects its two adjacent 
parabolas at (N + 1/2) $ (Fig- 0)- The current along the array is given by 

W = -^ = - r -(*otf-*erf). (88) 
It is proportional to the expectation value of the velocity of the charge soliton 

(*> = A^T = -W- (89) 



This is the quantum version of the relation (p2j) . We see that the external 
flux induces a persistent motion of the charge soliton, which is manifested in 
a persistent current along the array. As was shown above, no net number of 
flux quanta can tunnel in or out of the junction. However, during the motion 
of the soliton one can think of flux quanta flowing in and out of the array 
through the junctions, thus forming a flux loop around the moving center 
of the soliton. (A similar idea for 2-D superconducting films was given in 



32fl.) This interpretation is dual to the interpretation of the fluxon in a 



long Josephson junction as a (charge) current loop. The charge soliton's 
persistent current has the same origin as the persistent current of electron 
in a metal ring ||33|| . It is a manifestation of the Aharonov-Bohm effect [3| 



of a charged particle encircling a flux tube, and its persistency is due to the 
particle being in an exact eigenstate of the system. However, in contrast 
to the electron, the charge soliton is a macroscopic particle (Ac ~ 100 /im), 
so the possibility that it exhibits quantum effects is very intriguing. The 
quantum behaviour of the charge soliton is dual to the quantum behaviour of 



the fluxon in a long Josephson junction [19|. The latter exhibits a persistent 



motion in response to an external bias charge, which is the manifestation 



of the Aharonov-Casher effect [35|. Being a magnetic particle, this motion 
results in a persistent voltage across the junction. 

A weak spatial inhomogeneity in the array, e.g. non-identical grains or 
junctions or disordered grains, gives an additional X - dependent term in 



the Hamiltonian (82). The momentum is not conserved anymore, and flux 
quanta can tunnel across the array, reflecting in the spectrum by gaps which 
are opened at the intersection points of the parabolas (see Fig. [7]). If the 
array is now adiabatically biased by a time varied flux source, the persistent 
current oscillates as a function of $ ex .£ with a period $o- in each period a 
flux quantum tunnels across the array. This tunneling creates a current in 
the inverse direction to the existing current, thus eliminating the net current 
and reducing the energy. Since the energy bands are exact eigenstates, the 
tunneling process is a coherent one. When the external flux is not equal to 
an integral number of flux quanta, the quantum state of the array is a super- 
position of two flux quantum states. The amplitude of a persistent current of 
one charge soliton decreases as the amount of inhomogeneity increases. The 
maximal amplitude, corresponding to a vanishing amount of inhomogeneity, 
is of the order of 0.1 nA. 
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Figure 7: The spectrum of a quantum charge soliton in a 1-D ring-shaped 
array of serially coupled Josephson junctions as a function of an external flux, 
&ext- In an ideal ring the spectrum consists of inductive energy parabolas 
without a possibility of crossing at the intersection points. When there is 
some inhomogeneity in the ring (e.g. due to disorder), gaps are open at the 
intersection points, and the spectrum develops into energy bands. 

5.3 Other Quantum Effects 

The quantum nature of the charge soliton can be revealed in transport phe- 
nomena as well. For instance, if solitons are sent through a ring-shaped array 
connected to two leads (all consist of serially coupled Josephson junction), 
and the dephasing mechanisms are suppressed, we expect that they will split 
into partial waves propagating along the two arms of the ring. The partial 
waves will then interfere at the outgoing leads, with the interference pattern 
being dependent on the length of the arms and on an external flux applied 
through the center of the ring. The transmission of quantum charge solitons 
through the ring is thus expected to show oscillations as a function of the 
external flux and of the optical path similar to the h/e oscillations in the 
transmission of electrons through a metal ring fl36| , and in analogy to the 
transmission of fluxons through a Josephson junction ring JT7| . 

5.4 Dephasing Mechanisms 

The quantum phenomena described above were a consequence of the fact 
that in our approximation the Hamiltonian fl82|) was a one particle Hamil- 
tonian. Thus, even in the presence of a weak inhomogeneity, the degree of 
freedom associated with the charge soliton's center of mass (Xq) can main- 
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tain its quantum coherence. In order to make the model more realistic, one 
should take into account interactions between the soliton and other degrees of 
freedom. These interactions can produce, in principle, phase breaking mech- 
anisms. Whenever the phase breaking length, defined as the length over 
which the soliton's phase has an uncertainty of 2ir, is smaller than the length 
of the array, the quantum phenomena exhibited by the charge soliton will be 
suppressed. As in the case of the fluxon in a long Josephson junction |3Cj, we 
can distinguish between internal and external dephasing mechanisms. The 
internal mechanism is due to the interaction between the charge soliton and 
the other degrees of freedom of the junction, i.e. the plasmons. When the 
sine-Gordon model is exact and continuous, the system is completely inte- 
grable and the soliton is decoupled from the plasmons. Nevertheless, it has 



been shown in the context of the fluxon in a long Josephson junction |pC 
that there is a possibility of dephasing in this case as well. In order to avoid 
this dephasing, the temperature should be below the plasmons' energy gap. 
In the context of the charge soliton, where the sine-Gordon model is only 
an approximation and the system is discrete, we expect that the plasmons 
give rise to a stronger dephasing due to their inelastic interaction with the 
soliton. From the study of the discrete sine-Gordon model it is known that 
the rest energy of a soliton whose center resides in a junction is higher than 
the rest energy for a soliton whose center resides in the middle of a grain 
f39fl . Thus the soliton propagates in a periodic potential and not in a flat 
one. This deviation from the continuum model produces a coupling between 



the plasmons and the soliton. The soliton can emit or absorb plasmons |38 
|3Tfl , and the circulating soliton can become phase locked with this plasmons 
p0| . This effect has been recently observed for the fluxon in the discrete 
long Josephson junction [fy]]. We expect that similar phenomena occur in 
the system we study here when the continuum condition ( |26| ) does not hold. 
Apart from producing a phase breaking length, these phenomena will affect 
the classical dynamics as well, for instance by creating resonances in the I-V 
characteristic. The influence of both the discreteness of the array and the 
deviation from the exact sine-Gordon model on the classical and quantum 
mechanical dynamics of the charge soliton should be studied further. 

The most important external dephasing mechanisms are due to interac- 
tion with quasi-particles, which was neglected in our model. Since the bulk 
superconductors energy gap, A, is typically of the same order or higher than 
the plasmons' energy gap, the condition needed to suppress the thermal acti- 
vation of the plasmons is sufficient to suppress the thermal activation of the 
quasi-particles. The effects of thermal quasi-particles will be studied else- 
where. The quasi-particles can destroy the quantum coherence of the array 
in another way, which is temperature independent. The complete spectrum 
of a single junction includes charging energy parabolas associated with quasi- 
particles as well, which are separated in the charge axis by e. Excitation of 
quasi-particles leads to transitions between these parabolas, thus destroying 
the quantum coherence of q. This effect can be neglected if the charging en- 
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ergy of the quasi-particles plus the superconducting energy gap is larger than 
the Cooper pairs charging energy, i.e.: 2A + e 2 /8C > e 2 /2C or yA > E C - 
Since Ec should be smaller than A for the existence of the Josephson effect, 
this condition is met automatically. 

6 Summary 

We have studied a 1-D array of serially coupled Josephson junctions in the 
limit when the kinetic inductance of the superconducting grains dominates 
over the Josephson inductance. In this case the array is described by variables 
which are defined on the junctions and not on the grains. We have shown 
that the large kinetic inductance decouples the junctions quantum mechani- 
cally. As a result each junction is characterised by a periodic charging energy. 
This periodic energy, when combined with the inductive energy of the grains 
and the charging energy between the grains and the substrate, gives rise to a 
model with topological solitons excitations. Thus we have found that an ex- 
cess Cooper pair in the array creates a charge soliton via polarization of the 
superconducting grains. The charge soliton is a dual topological excitation 
to the fluxon in a long Josephson junction. We have studied the classical 
dynamics of the charge soliton, and shown that in the presence of dissipation 
and an external voltage the I-V characteristic of the array should consist of 
saturation branches corresponding to the number of charge solitons in the ar- 
ray. We have quantized the charge soliton semi-classically, showing that this 
quantization is consistent with the large kinetic inductance. We have found 
that a quantum soliton in a flux-biased ring-shaped array is expected to show 
persistent motion, manifested in a persistent current. A weak inhomogeneity 
in the array gives rise to a coherent current oscillations. These phenomena, 
which are usually associated with electrons (or Cooper pairs) suggests that 
the quantum charge soliton can be considered as a macroscopical quantum 
object. Finally, we have discussed possible internal and external dephasing 
mechanisms of the charge soliton. These mechanisms deserve a future study. 
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